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Abstract 

We use a discrete worldline representation in order to study the con¬ 
tinuum limit of the one-loop expectation value of dimension two and four 
local operators in a background field. We illustrate this technique in the 
case of a scalar field coupled to a non-Abelian background gauge field. The 
first two coefficients of the expansion in powers of the lattice spacing can 
be expressed as sums over random walks on a d-dimensional cubic lattice. 
Using combinatorial identities for the distribution of the areas of closed 
random walks on a lattice, these coefficients can be turned into simple 
integrals. Our results are valid for an anisotropic lattice, with arbitrary 
lattice spacings in each direction. 


1 Introduction 

1.1 Motivation 

The classical statistical approximation (CSA) is an approximate scheme to study 
in real time the dynamics of a system of helds, as an initial value problem. It 
has been used in cosmology [1, 2, 3], in high energy nuclear physics to study 
heavy ion collisions [4, 5, 6], in studies of the Schwinger mechanism [7, 8], and in 
cold atom physics [9, 10]. In this approximation, the time evolution of the fields 
is classical (i.e. deterministic), and one averages over fluctuations of their initial 
conditions. Obviously, this scheme neglects all the quantum effects that would 
normally affect the time evolution of a system. In the path integral language, 
it corresponds to taking the saddle point of the integral. This can be justified 
if the typical field amplitude in the system under consideration is large, so that 
the commutator of a pair of fields is much smaller than the typical held squared. 
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Some quantum effects can nevertheless be included in the classical statistical 
approximation via the fluctuations of the initial fields. Indeed, it can be shown 
on general grounds that the leading (i.e. of order h) quantum effects come 
entirely from the initial condition for the density operator of the system, while 
the quantum corrections that alter its time evolution only start at the order 

In fact, there is a unique statistical ensemble of initial classical fields such 
that the CSA coincides with the exact 0{h) result for all observables [1, 11, 12]. 
These initial fields have a flat spectrum in momentum space, that extends to 
arbitrarily large momenta. 

Because the CSA implemented in this manner^ contains all the 0(h) con¬ 
tributions of the underlying theory, it also contains all their singularities, and 
in particular the ultraviolet divergences. In addition, it contains some of the 
higher order contributions, but not all of them since the quantum corrections 
to the time evolution are missing in this approximation. It has been shown 
recently that this leads to a dependence on the ultraviolet cutoff that cannot 
be disposed of by the usual renormalization of the parameters of the theory. 
This can be seen in computer simulations using the CSA [15], by a perturbative 
analysis of the graphs that arise in the CSA [16], and also from the study of 
cutoff effects in the classical approximation of the Boltzmann equation [17]. 

At the moment, it is not known whether the CSA can be modified in order 
to remove these unwanted terms. However, regardless of this interesting theo¬ 
retical question, it is important to have a good understanding of the structure 
of the standard 1-loop ultraviolet divergences. Indeed, since they are identical 
in the CSA and in the exact theory, they can be removed by the usual renor¬ 
malization procedure. But their form may be quite complicated in the lattice 
formulation of the CSA, especially for a generic lattice that may have anisotropic 
lattice spacings^. Generically, this requires the calculation at one-loop of the 
expectation value of interest, with lattice regularization and propagators, in the 
presence of a non-Abelian background field. Unfortunately, lattice perturbation 
theory is quite complicated, even for this seemingly simple task (see [18] for a 
review). The main issue is the treatment of the background field, and the fact 
that one recovers gauge invariant results by combining several pieces that are 
not individually gauge invariant (see the appendix E for an example of such 
calculation in the present context). 

In the present paper, we pursue a different approach in order to obtain 
these 1-loop quantities, based on the so-called worldline formalism. Historically, 
this formalism emerged from ideas based on the limit of infinite tension in 

^Another common implementation of the CSA is to use a spectrum of field fluctuations 
with a compact (or at least falling faster than 1/fc) momentum spectrum [4, 5]. This type of 
initial distribution corresponds to a classical ensemble of quasiparticle excitations, instead of 
quantum fluctuations. This version of the CSA is free of any ultraviolet divergences (for a 
spectrum that falls like k~^, one gets some ultraviolet divergences, but the resulting approx¬ 
imation is super-renormalizable [13, 14]), but it also does not contain anything quantum. It 
may coincide with the underlying theory at tree level, but not beyond. 

^In applications to heavy ion collisions, it is common to have a much smaller lattice spacing 
in the direction of the collision axis. 
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string theory^ [20, 21], and it soon appeared that it provides a powerful way 
of organizing field theory calculations, especially when gauge symmetry is part 
of the game. A pure field theory understanding of this formalism was later 
proposed in ref. [22], by a method which is similar to Schwinger’s proper time 
representation. For reviews on this approach, the reader may consult refs. [23, 
24]. This formalism has been applied to the evaluation of effective actions 
[25, 26, 27, 28], to the study of pair production by an external field [29, 30] or 
the Casimir effect [31, 32]. Numerical algorithms based on this formalism have 
also been proposed [33, 34, 32]. Since our goal is to apply this formalism to a 
lattice field theory, the closest work we are aware of is in refs. [35, 36], where a 
new method for evaluating functional determinants in terms of worldlines was 
proposed. 

In the present paper, we use this formalism in order to obtain useful expres¬ 
sions for 1-loop expectations values in a lattice field theory coupled to a (fixed) 
gauge background. As we shall see, the worldline formalism is well suited for 
this application because it enables one to have only gauge invariant objects at 
all stages of the calculation. Then, we use these expressions in order to study 
the limit of small lattice spacing. In this limit, we obtain an expansion in terms 
of the background field strength, with coefficients that are given by sums over 
closed loops on the lattice, weighted by powers of the area enclosed by the 
loop. Thus, the worldline formalism relates the coefficient of the short distance 
expansion to combinatorial properties of closed loops on a cubic lattice. 

In the simple case where the lattice spacings are the same in all the di¬ 
rections, the combinatorial formulas we need were already known and can be 
found in ref. [37]. The generalization to anisotropic lattice spacings requires 
some combinatorial formulas that we could not find in the literature. A numer¬ 
ical exploration of all the random walks of length < 20 led us to conjecture a 
number of such formulas (discussed in the appendix C), that extend and gen¬ 
eralize those of ref. [37]. A proof of these formulas^ is presented in a separate 
paper, ref. [38]. 

1.2 Model 

In order to keep things rather simple and focus on the main aspects of the 
worldline formalism, we consider a complex scalar field coupled to an external 
non-Abelian gauge field. This background field is given once for all, and does not 
fluctuate. We neglect the self-interactions of the scalar field. The Lagrangian 

® An earlier example of this approach is a string-inspired calculation of the 1-loop /3 function 
of Yang-Mills theory [19]. 

■^Eqs. (142) and (146) can be viewed as “unsummed” versions of the formulas (1.5) and 
(1.6) of ref. [37], while eqs. (143), (144), (145), (147), (148) and (149) seem to be totally new. 
The appendix E provides an indirect proof of eq. (142), since we rederive the expansion of 
(0*(0)</)a(0)) -that relies on eq. (142) in the worldline approach- using lattice perturbation 
theory. 
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reads : 


d 

= (/?''</>), ( 1 ) 
M=1 

where — igA^^ is the covariant derivative in the presence of the back¬ 

ground field. We also assume that the system is initialized at = —oo into 
the perturbative vacuum. 

We consider the expectation values of local gauge invariant operators made 
of the field </>, e.g. 0*^, {D^(j))*{Dy(j)). With the Lagrangian given 

in eq. (1) and the vacuum state as initial condition, these expectation values 
are given by a 1-loop graph in a background field. These loops contain a pure 
vacuum contribution which is ultraviolet divergent. In addition, for operators 
that have a sufficiently high dimension, there may be subleading ultraviolet di¬ 
vergences, whose structure is however strongly constrained by gauge invariance. 
Our goal in this paper is to investigate the structure of these divergences, with 
a lattice regularization. 

We work with an Euclidean metric, in d space-time dimensions, on a discrete 
cubic lattice. For definiteness, the spatial directions are chosen to be 1, • • • , d — 
1 and the direction d can be considered as the time direction (although this 
distinction is hardly relevant with an Euclidean metric). Our goal is to study 
the general case of arbitrary lattice spacings ai, • • • , in each direction, but 
in the sections 2 and 3 we expose the formalism with an isotropic lattice for 
simplicity. The coordinates are labeled xi to Xd-, and we denote by 1, • • • ,d the 
vectors corresponding to one lattice spacing in each of the directions. Therefore, 
any point on the lattice can be represented as x = a;il-|-----l- where the 
Xi are integers. 


1.3 Preview of the results 

In d dimensions and for completely arbitrary lattice spacings in each direction, 
the expectation value of |0(O)p at I-loop in a background field admits the 
following expansion in powers of the lattice spacings, 


(C(O)0a(O)) ^ = 


{a.^0} 2dntia* 


t^adj (1) ^0 




i<j 


with 


and 


Cn = 


/^ij;ij _ hihj 


dte-^ 

poo 




( 2 ) 

(3) 

(4) 
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Iq and Ii are modified Bessel functions of the first kind. In all these equations, 
we denote 


a = - 


1 

1 — 1 



(5) 


These formulas are the archetype of the results obtained in this paper. We also 
derive similar formulas for bilocal operators of the form (^(j)*g^{Q)'Wab{lxo)4’b{x)), 
where the point x is separated from the origin by 1 or 2 lattice spacings. In 
these operators, Wat( 7 ^ 0 ) is a Wilson line along a path 73,0 connecting x to 0, 
which is needed in order to have a gauge invariant operator. We shall see that 
the leading term of the expansion in powers of the lattice spacings does not 
depend on the choice of this path, while the second term in general depends on 
this choice. 

In the rest of this paper, we use the lattice worldline formalism in order 
to demonstrate these formulas for the coefficients of the expansion. We hrst 
obtain intermediate representations of these coefficients in terms of sums over 
all the closed random walks on the lattice, which relate their values to some 
combinatorial properties of random walks. These formulas can then be trans¬ 
formed into the integral representations listed above, by using the 2 -dimensional 
combinatorial formulas of the appendix C. 


1.4 Outline of the paper 

In the section 2, we derive in detail the worldline formulation of the expectation 
value ((/)*(O) 0 a(O)) on a lattice with isotropic spacings, and its short distance 
expansion. A subsection is devoted to the discussion of infrared divergences and 
their manifestation in the worldline formalism. We also introduce in this section 
the Borel transformation that turns the combinatorial sums into integrals. In 
the section 3, we extend this study to bilocal operators, i.e. operators that 
contain a <f)* and a (j) evaluated at separate lattice spacings. This extension is of 
great practical importance, because these operators appear in the discretization 
of covariant derivatives. The section 4 generalizes all the previous results to 
a more general lattice setup, where each direction of space-time has its own 
lattice spacing. As an illustration, we study the limit where one of the lattice 
spacings is much smaller than the others, and we apply this to a discussion of 
the energy-momentum tensor. The section 5 is devoted to concluding remarks. 

A number of more technical aspects are discussed in several appendices. In 
the appendix A, we show how this formalism is modified on a finite lattice with 
periodic boundary conditions (in the main part of the paper, we take the limit 
of zero lattice spacing at fixed physical volume, so that the size of the lattice 
becomes infinite and the boundary conditions are irrelevant). In the appendix 
B, we derive the leading coefficient of the expansion for bilocal operators with 
an arbitrary separation between the two fields. The appendix C discusses all the 
combinatorial formulas that are necessary in the case of anisotropic lattices, and 
in D we recall the connection between the statistics of the areas of closed loops 
on a two-dimensional lattice and the spectral properties of the so-called almost- 
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Mathieu operator. In the appendix E, we obtain the short distance expansion of 
((/)* (0)(/)a(0)) from lattice perturbation theory, mainly to illustrate the technical 
complexity of this approach. In the appendix F, we study from the outset a 
hybrid description in which one of the directions (e.g. time) is treated as a 
continuous variable, while the others remain discretized and we show that this 
is equivalent to starting from a fully discrete description and taking one lattice 
spacing to zero. 


2 Local operator <|^0*(O)0a(O)) 

2.1 Discrete heat kernel 


In order to establish the formalism, consider first the expectation value of the 
operator (())*(0)(/)a(0)) (the color indices are summed to ensure gauge invariance) 
at one loop. To lighten the notations, we evaluate the expectation value at the 
point = 0, but all our considerations are completely general. Since this 
expectation value is given by a 1-loop graph, we can first write 

(</>a(0)(('a(0)) = {x>^ = 0\-^\x>^ = 0) . (6) 


The standard heat kernel approach would be to write 


1 


ds exp{—sD^) . 


(7) 


However, in our case it is more convenient to use a discrete version of this 
formula^, 

^ = f;(i-a»Dy2dr. (s) 

n—0 

which is also exact. We have multiplied by the lattice spacing squared in 
order to get a dimensionless combination. The purpose of the factor 2d (where 
d is the number of space-time dimensions) will become clearer later on. 

Consider now a sequence of functions Pn(x) defined on the lattice, and sat¬ 
isfying the following iteration rule: 


Pn+i = (1 - a'^D'^/2d) Pn . 


(9) 


If we define Po(x) = Sx, then we have 

(0|(1 - a^D‘^/2d)^\0) = P„(x) . (10) 

If we interpret Poix) as a probability distribution localized at the point x^ = 0, 
then P„ is the probability distribution after n iterations of the process described 
in eq. (9). In other words, it is the probability that this process starts and returns 
at the point x^ = 0 after exactly n steps. 

®Up to a rescaling, the integer n is a discrete version of the integration variable s. 
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2.2 Vacuum case 

Eq. (9) may be rewritten as 


(,2 

P„+i-P„ = -—(11) 

If we view the index n as a discrete fictitious time, and if the metric is Euclidean, 
then this is a discrete diffusion equation and the evolution of the probability 
distribution P„ can be remapped in terms of random walks. 

Eor illustration purposes, consider first the free case. We have 





-/(*,•••) 


/(z + 1, ■ ■ ■) + /(f — 1, ■ ■ ■) 

2d 


( 12 ) 


The eq. (11) can then be written more explicitly as 


Pji+l i'i: ‘ ‘ ‘ 


Pn(i + 1, • • • ) T Pni'i ~ 1) • ■ ■ ) 
2d 


(13) 


where the sum in the numerator extends to all the nearest neighbors. This 
equation describes a random walk where at each step one moves to one of the 
adjacent sites of the lattice with probability l/2d. We see now the reason for 
the peculiar normalization® in eq. (6) : by doing this, we can eliminate the 
term proportional to P„(z, • • •) in the right hand side, i.e. the possibility for the 
random walk process to stall during the step^. 

P„(0) is the probability that such a random walk returns at the point 0 
after n steps. Geometrically, this means that the random walk is a closed loop 
of length rz. Since at each step, there are two possibilities to move in each 
direction, the total number of random walks of length n is (2d)". Pn{x) is thus 
the number of closed random walks of length rz, divided by the total number 
(2d)". Therefore, we can write 


(C(0)((-a(0)) = Y. (1) > (14) 

"=0 '' ’ 76 r 2 „( 0 , 0 ) 

where r 2 „( 0 ,0) is the set of all closed random walks of base point 0 (i.e. starting 
and ending at 0) and length 2n on the lattice (the length of such a closed path 
must be an even number). A few of the closed paths involved in eq. (14) are 
illustrated in the figure 1. 

In the vacuum, the double sum is independent of the lattice spacing. It is 
just a pure number that sets the normalization of the result. The trace in the 
adjoint representation comes from the summation over the color indices, and 
brings a factor — 1. 

® Alternatively, one could view this normalization as choosing a specific ratio between the 
size of the steps in the fictitious time and the lattice spacing. 

^By excluding the possibility that the random walk stalls, we ensure that the number of 
steps n is also the length of the path. 
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Figure 1: Example of closed random paths on a 3-dimensional cubic lattice. 
The blob indicates the location of the base point 0. 

2.3 Non-zero background field 

On the lattice, the background field is represented in terms of compact link 
variables Ui(x) in order to preserve an exact gauge invariance despite the dis¬ 
cretization. In terms of these link variables, the covariant derivative squared 
becomes. 


••) = -/(*,•••) 

+■ - + • • • . (15) 

Therefore, when the links are not unity, the random walk is biased by the 
background field. The end result is that eq. (14) is modihed into 

1 °° 1 

"=0 ^ '' 76r2„(0,0) 

In words, the SU{Nc) identity matrix in eq. (14) is replaced by a Wilson loop 
14 ^( 7 ) obtained by multiplying all the link variables along the closed contour 
7 . This formula is manifestly gauge invariant, since Wilson loops are gauge 
invariant. Note also that this formula is exact at 1-loop on the lattice. 

2.4 Notations and basic facts about closed random walks 

In the previous subsections, we have introduced r2„(0,0), the set of all the paths 
of length 2n drawn on the lattice, with endpoints 0 and 0 (i.e. closed paths). 
More generally, we will denote r„( 0 ,a:) the set of paths of length n from 0 to 
a point X (all these paths have the same parity, which is also the parity of the 





































sum of the coordinates of the point x). To avoid encumbering the notation, 
we do not specify the dimension of the lattice in which these paths should be 
considered, since the context of the formula in which the notation appears is 
sufficient to make this obvious. 

When we need an explicit representation for a path, we denote it by the 
sequence of the hops it contains, such as 

7 = (17) 

(read from left to right.) The notation 3“^ denotes a hop in the —X 3 direction. 
The empty path will be denoted 7 = 1 , and the concatenation of two paths 71 
and 72 is denoted by 7 = 71 ( 8)72 (read again from left to right, so that 71 is the 
first part of the resulting path). Obviously, 107 = 7(8)1= 7. 

In two dimensions, we will also introduce later in the paper the subset 
r„i_„ 2 ( 0 ) 2 ;), made of all the paths connecting 0 to x and making exactly ni 
hops in the +xi direction and n-i hops in the +X 2 direction. The numbers n'l 2 
of hops in the opposite directions, —Xi and —X 2 , do not need to be specified 
explicitly since it can be inferred from ni _2 and the coordinates of the point x. 
Indeed, if x = xil + X 22 , we have 


n\ = ni — Xi , n '2 = n 2 — X 2 ■ 

Obviously, r„( 0 ,x) and ( 0 : 2 ;) are related by 


r„(o,x) 


Tni ,712 (®i 2)) . 


ni+n2 = 


n-\-xi+X2 

2 


(18) 


(19) 


In order to develop some intuition with formulas such as eq. (16), let us 
recall here some elementary properties of closed random walks. Let us consider 
a closed random walk made of 2 n hops, all illustrated in the figure 2 in two 
dimensions. For such a random walk, one has the following properties: 

i. the length of the path is obviously 2 na, 

ii. the typical size of the domain explored by the random walk grows only 
like ^/na, 

iii. the area enclosed by the random walk (or the area of its projection on a 
plane in d > 2 dimensions) grows as na^. 

The property ii plays a role in the infrared behavior of the quantity under 
consideration. Indeed, as we shall see later, infrared singularities arise when the 
contribution of “large” random walks does not decrease fast enough. Similarly, 
iii plays a role in the second term in the expansion in powers of the lattice 
spacing. 
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Figure 2: Closed random walk on a square lattice in two dimensions. The length 
of the path is 2 na, while the diameter of the domain explored by the random 
walk is of order y/na and its area is of order na^. 

2.5 Continuum limit a —)■ 0 

If we let the lattice spacing a go to zero, while the background is held fixed in 
physical units, the Wilson loop ^(y) goes to the identity because the closed 
loop 7 shrinks to a tiny loop of base point 0. Therefore, it can be approximated 
by the exponential of the magnetic flux across a surface E whose boundary is 7 
(this surface is a tiling of elementary lattice squares), 

W(7) exp [tga^ ^ A^,{^)Fr{0r} , (20) 

where is the algebraic area, measured as a number of plaquettes since we 

have already pulled out a factor a^, of the domain enclosed by the projection of 
the contour 7 on the (/r, v) plane. The orientation of 7 dictates the orientation 
of the projection, which in turn controls the sign of ^^^( 7 ). Note that in this 
limit, the field strength can be considered uniform across the entire lattice, and 
therefore is evaluated at the point = 0 . 

Since the loop size tends to zero when a —>■ 0, we can do a Taylor expansion 
of the exponential. In order to get a non-trivial answer after taking the trace, 
we must go to second order: 


tradj (W(7)) « tradj ( 1 ) 


(A factor 1/2 comes from 


p<(7 

( 21 ) 

tradj(i“i^) = By plugging this in the formula 
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(16), we obtain the following expansion 


(C(0)<(<a(0)) 


1 

2 d 


fl'adj (l) 




4 


^Fr(o)i^r(o)c" 

p<<7 


(22) 


where the coefficients Co and are purely geometrical quantities defined 

by sums over all the closed loops on the lattice 


^ (2d)2" ^ ^ 

n -0 ^ ^ 7er2^(0,0) 

OO -J 

^4 = X/ f2d')2n X/ ■ (23) 

"=0 ^ ' 76r2„(0,0) 


(Note that in these formulas, we have replaced n —> 2n since only random 
paths of even length can be closed.) The equation (22) provides an explicitly 
gauge invariant expansion in powers of the lattice spacing. The coefficients 
(that remain to be calculated) are geometrical quantities that depend on the 
dimension and the lattice under consideration, but not on the background field. 

In the second of eqs. (23), Ap^{j) is the area of the surface enclosed by 7 
projected on the pv plane. Several remarks are in order about this quantity: 

• These areas are “algebraic”, in the sense that they may have a sign that 
takes into account the orientation of the boundary, and a multiplicity that 
depends on the winding number. 


• There are many surfaces with the same boundary 7 . does not 

depend on this choice but only on the boundary. 


• They do not depend on the base point = 0. Specifying a base point is 
only necessary in order to avoid counting multiple times loops that have 
the same shape up to a translation. 


2.6 Zeroth order coefficient 


It is possible to provide an integral expression for the coefficient Cq, starting 
from the combinatorial formula that explicitly counts the number of closed 
random walks in terms of the number of hops in the d directions (respectively 
2ni,2n2,--- ,2nd), 


f. ( 2 n)! 


E 

iilH-|-nd=™ 


1 

m!^ • • • nd!^ 


(24) 


The factor (2n)! prevents the complete separation of the sums over the n^. 
However, it can be removed by a Borel transformation : 

cOO 

Co= dte-^Adii^), (25) 

do 


II 



where we denote 


A,{x) ^ 5] 


n=0 niH- \-nd—n 

-\ d 

d 


ni!^ • • • nd!^ 


°° a;2p 


E X 

pI2 


.p=0 


= m ^), 


(26) 


where Iq is a modihed Bessel function of the hrst kind. Therefore, we have 


Co = 


nOO 

^0 


In 3 dimensions, this leads to an explicit formula ([39]-§6.612.6) : 

v^- 1 


Co = 


« 1.51638606, 


d=z 327r3 

while in 4 dimensions we have only been able to evaluate it numerically. 


(27) 


(28) 


Co « 1.23946712 . 

d=A 


(29) 


2.7 Variance of the areas of closed random walks 

In the second term of the expansion in powers of the lattice spacing, we need 
the quantity 


OO ^ 

C 4 ^ ^ ■ (30) 

"=0 i i 76r2„(0,0) 

A central result for the rest of our discussion in the isotropic case is the value of 
this sum in two dimensions. In d = 2, the variance of the algebraic areas enclosed 
by closed random walks of length 2n is given by the following formula [37] 
(eqs. (1.4)-(1.5)) : 


E (^12(7))^ 

. 76 r 2 „( 0 . 0 ) 


dim 2 


/2 n\ ^ n^(n — 1 ) 
\ n / 6 ( 2 n — 1 ) 


(31) 


This formula is all we need in order to evaluate the coefficient C 4 for d = 2. But 
note that this coefhcient diverges in d = 2: using Stirling’s asymptotic formula 
for the factorial, one can see that the sum over the length 2 n of the path is 
divergent (in this formalism, this is the counterpart of an infrared divergent 
loop integral in low dimension). 

In higher dimensions, the hrst thing to notice is that it is sufficient to consider 
jjL = p,v = a (if there is a mismatch of the indices, the average over all closed 
loops gives zero since the area is signed). For the sake of dehniteness, let us 
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Figure 3: Example of closed random walk on a 3-dimensional cubic lattice, with 
2n = 8 hops in the 12 plane and 2m = 4 hops in the 3-rd direction. The 12 
plane is highlighted in orange, and the projection of the closed loop on the plane 
is shown in gray (it has an area A 12 = —2 in this example). The blob indicates 
the location of the base point 0 . 


choose fj, — 1, ly = 2. For a given closed loop 7, the area ^12(7) is the area of its 
projection on the 12 plane. Every closed random walk in d dimensions can be 
decomposed into hops that are in the 12 plane, and hops orthogonal to this plane 
(see the figure 3). The latter disappear in the projection on the 12 plane, and 
therefore do not play any role in the calculation of the area ^ 12 ( 7 ). Moreover, 
the projection of 7 in the 12 plane is itself a closed loop in 2 dimensions, while 
the sequence of the transverse hops is a closed loop in d — 2 dimensions. Let us 
denote 2 n the number of hops in the 12 plane and 2 m the number of hops in 
the transverse directions. One can rewrite the coefficient C 4 as follows 

^I 2 ;i 2 _ 1 /2(m-|-n)\ 

4 - (2d)2("+™) V 2m / 


X 

E 1 

X 

E (^ 12 ( 7 ))' 


_<TGr2„(o,o)_ 

dim d-2 

7Gr2„(0,0) 


The binomial factor in the first line counts the number of ways to intertwine 
the 2 m transverse hops and the 2 n in-plane hops. In the second line, the first 
factor is the number of length 2 m closed random walks in d — 2 dimensions, and 
the second factor is the squared area summed over all closed loops of length 2 n 
in two dimensions. This latter factor is given by eq. (31). For the first factor, 
since we are interested primarily in d = 3 (d — 2 = 1) and d = 4 (d — 2 = 2), we 
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can use the following standard results 


E 1 

crer 2 ™(o,o) 


crer 


E 1 

(0.0) 


J dim 1 


/ 2 m 
\ m 


2 m 


(33) 


J dim 2 


Therefore, for these dimensions, we find the following expressions for the coef¬ 
ficient 


= V 


d=Z 


(=0 


(2/)! (2n)! n^(n—1) 

6 ^^ {I — n)!^n!4 6(2n — 1) ’ 


and 


C 


I 2 ;i 2 _ (20- ~ n))!(2n)! n^(n — 1) 

d=4 8^* il — n)\^n\^ 6(2n — 1) 

1=0 n=o ' 


(34) 


(35) 


For a given total length 21 = 2m + 2n of the random walks, the variance of 
the projected algebraic area is thus expressed as a sum of I -I- 1 terms, whose 
evaluation is very easy (especially compared to a direct evaluation by exhausting 
the list of all random walks of this length, since there are {23)^^ such walks). In 
the table 1, we list as an example the summands (at fixed 1) for the 3 dimensional 
case (eq. (34)). We have performed an exhaustive search of all the closed random 
walks up to 21 = 14, and we have checked the agreement between this direct 
computation and the formula (34). The values listed for 21 > 14 were solely 
obtained from eq. (34). 

It is also possible to establish an integral representation of valid in 

any dimension d, similar to eq. (27) for Cq- The first step is to introduce the 
combinatorial representation for the factor that counts the closed random walks 
of length 2m in d — 2 dimensions. This leads to 


OO 


1 (2(m-|-n))! (2 to)! 

(2(i)2(™+”) (2m)!(2n)! ^ mi\'^ ■ ■ ■ md- 2 ' 

m,n —0 ^ ^ \ / \ / mi-\ -hmd_2=m 


12 


(2n)!^ n?{n — 1) 


n!"* 6(2n — 1) 

We can separate the sums over m and n by a Borel transformation, 


^ 12;12 _ 
— 




(36) 


(37) 


with 

Cd{x) = 


“ ^ 2 (m+n) 

(2m)!(2n)! ^ mi!^ • • • md-2!^ 

m,n—0 ^ / \ / miH- \-md- 2 —^ 


- It\2x) lH2x) . 


(2m)! (2n)!^ n^(n — 1) 

n!^ 6(2n — 1) 

(38) 
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21 

# (paths) 

# (closed paths) 

(0,0) (^12(7)) 

q 21 

2 

36 

6 

0.0000000000 

4 

1296 

90 

0.0061728395 

6 

46656 

1860 

0.0102880658 

8 

1679616 

44730 

0.0133363816 

10 

60466176 

1172556 

0.0158369532 

12 

2176782336 

32496156 

0.0180064306 

14 

78364164096 

936369720 

0.0199490385 

20 

3fP 


0.0249038527 

100 

36®° 


0.0600254031 

200 

36100 


0.0856471034 

1000 

36500 


0.1928668060 

2000 

361000 


0.2729940025 


Table 1: Exhaustive enumeration of random walks and closed random walks in 
3 dimensions, up to the length 21 = 14. The last column gives the corresponding 
contribution to . The values for 21 > 14 are obtained from eq. (34). 


Therefore, we have 


^ 12;12 

'-"4 


1 



dte-^ 


(39) 


2.8 Infrared divergences 

As one can see in the table 1, the summands do not decrease at large path 
lengths I, and the sum over I is divergent. In d = 3, the summand grows as 
for large path lengths, and it goes to a constant in d = 4. If we cut off the sum 
over I at some ?max, this implies that 


^ 12;12 . 3/2 ^ 12;12 . 

'-^4 *max J *-^4 *max • 

a—3 a—4 


(40) 


This divergence for random walks that explore large regions in spacetime is 
the manifestation in the worldline formalism of an infrared singularity. This is 
corroborated by the fact the divergence is milder in d = 4 compared to d = 3. 
In the integral representation (39), this singularity appears as a divergence of 
the integral at large t : using the fact that /„(t) ~ at large t, we see that 

the exponential factors cancel and that the remaining algebraic factors decrease 
fast enough for convergence only if d > 6. 

In order to further investigate this, let us add a mass term® to the Lagrangian 

®One may view this mass as a temporary regulator for the infrared sector. Note that if 
we do not expand the observable in powers of the background field, then the background 
field itself would provide a natural infrared cutoff. The sum over I would naturally be cutoff 
when random walks reach a size comparable to the coherence length of the background field. 
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of the scalar field, 


d 

L = Y.^D^4’r{D^cf)-m^cj,*(j>. (41) 

^l=l 

In order to arrive again at a sum of random walks without stalls, one should 
start from the following formula 


2d 

a?(D‘^ + w?) 


Ea 


„2 D'^+rrd \n 
“ 2d > ’ 


n—0 


(42) 


where we have defined d = d+ \rn?a?. Most of the discussion is unchanged and 
eq. (16) becomes 


(<^:(o)</),(o)) 


1 

2Ja^-2 


OO 


E 

n—O 


1 

(2(i> 


tradj (W(7)) • 

76 r„( 0 , 0 ) 


(43) 


A crucial difference is that each hop in the random walk is now weighted by a 
factor l/2d instead of l/2d. Since d > d, this leads to an exponential reduction 
of the contribution of long random walks. The eqs. (34) and (35) are modified 
into 


OO 


/-rl2;12 

^4 


d=3 


E 

1=0 


{2iy- 

(6 + 


i 


E 

n—O 


(2n)! n^(n—1) 

{I — n)Pn\‘^ 6(2n — 1) ’ 


(44) 


and 


I 2 ;i 2 _ (20- {2(1 — n))\{2n)l {n — 1) 

(i= 4 (8 + 2^ (I — n)\^nl^ 6(2n — 1) 


(45) 


These discrete sums also have an integral representation in terms of modified 
Bessel functions, 


^ 12;12 

*-"4 


1 

1 ^ 





(46) 


The derivation of this formula follows the same line as that of eq. (39). 

The modification of the behavior for large random walks is readily seen. 


1 _ 1 1 
(2(1)2* “ (2d)2* + 


1 

(^ 




(47) 


The regularizing effect of the mass becomes effective for lengths I > d/m'^a^. 
This corresponds to random walks that explore a domain of size r > \fd/m, i.e. 
of the order of the Compton wavelength. The values of the summand in d = 3 
are shown for = 0.01 in the figure 4, and compared to the massless case. 
With a non-zero mass, the summand increases until it reaches a maximum and 

For instance, if the background field is incoherent beyond the length scale Q~^, then values 
I > (Qa)~^ are suppressed. 
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d = 3 



Figure 4: Summand in 3 dimensions, in the massless case (black curve) and for 
= 0.01 (red curve). 


then decreases exponentially, which ensures the convergence of the sum over 1. 
Note that in doing this, we do not expand the a dependence that comes from 
2d + m^a?. 

After this infrared regularization, the singular behaviors of eq. (40) would 
be replaced by 

C 4 - (ma) , C 4 ’ - [ma) . (48) 

a—3 a—A 

These formulas contain inverse powers of the lattice spacing. This tells us that 
the expansion of eq. ( 22 ), in which the second term is suppressed by a^, is upset 
by infrared singularities. In fact, the second term does not vanish when a 0, 
but instead is a term of order a^, both in d = 3 and d = 4, 

(</>a(0)</>a(0)) ^ = ^ . (49) 

The same is true when the infrared regularization is provided by the back¬ 
ground field rather than by a mass. In this case, the mass is replaced by the 
coherence scale Q of the background field in the above counting. Note that a 
proper calculation of the terms in requires a non-perturbative treatment of 
the background field (as we have seen before, it was the expansion in powers of 
the background field that caused the infrared divergence in the first place). 
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3 Bilocal operators (0*(O)Wa6(7xo)06(2^)) 


3.1 Worldline representation 

Up to now, we have considered only the local operator (f)* (0)(f)(0). However, 
since derivatives are represented on the lattice as finite differences, we need 
also to consider composite operators made of two elementary fields evaluated at 
different lattice sites. Let us illustrate this by the discretization of the operator 
(f)*D^D^(f>. On the lattice, this can be represented as follows, 




(50) 


Here, one of the derivatives has been discretized as a forward derivative and the 
other as a backward one. We have used covariant derivatives for the operator 
to be gauge invariant. After discretization, this leads to link variables connect¬ 
ing the two lattice sites where the scalar field and its complex conjugate are 
evaluated. 

All the terms in the right hand side of eq. (50) belong to a class of operators 
that contain two fields (f)* ■ ■ ■ (f) linked by a Wilson line, 

(</>:(0)Wafc(7xo)<(';.(a:)) • (51) 

In this equation, is a path (drawn on the edges of the lattice) connecting 
the points x and 0 where the two fields are evaluated. The Wilson line and the 
summation over the color indices ensure that this operator is gauge invariant. 
However, the choice of this path is arbitrary and is therefore a part of the 
definition of the operator under consideration. Except when the background 
field is a pure gauge, different paths correspond to operators that have distinct 
expectation values. 

It is quite straightforward to generalize the derivation done in the section 2 
in order to obtain the lattice worldline representation for this type of operator, 


{(l>*a(0)Wab(lxo)4>b(x)) 


1 


OO 


E 

n—Q 


1 

(2d)" 


^ tradj (W(7a;0(8)7)) 
7er„(0,a;) 


(52) 


In this formula, r„(0,x) denotes the set of all the random walks of length n 
that start at the point 0 and end at the point x, and 7a,o <8) 7 denotes the closed 
path obtained by concatenating one of these paths and the contour used in 
the definition of the operator, as illustrated in the figure 5. 


3.2 Continuum limit 

Let us now consider the continuum limit of these expectation values. There are 
two non-equivalent ways to view this limit: 
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Figure 5: Example of random paths that appear in the worldline representation 
of bilocal operators, on a 3-dimensional cubic lattice. The blobs indicate the 
locations of the base points 0 and x (here separated by 1 lattice spacing). The 
black link is the Wilson line W^'jxo) inserted between the two fields for gauge 
invariance. 


i. Keep the separation x-0 between the two points fixed in absolute units. 
When a —>■ 0, this interval becomes infinite in lattice units. This leads to 
milder ultraviolet divergences, since when the a —>■ 0 limit is performed in 
this way, one is in fact considering a non-local operator. 

ii. Keep x and 0 at fixed locations on the lattice while a —)■ 0. Therefore, the 
spacing a :-0 goes to zero in absolute units, and the limit o —)■ 0 corresponds 
to a local composite operator. 

The limit ii is the one we consider here, since we are interested in operators 
such as the one in eq. (50). The eq. (22) becomes 


{(j)l{o)yVabi.'ixo)Mx)) 


1 

2dad-^ 


-^0,7j,o f^adj (l) 

^ Fr{x)Fr{x) 

p<a 


,(53) 


where the coefficients £> 0 , 7^0 -^4 7 ^ 0 °^ generalizations of the coefficients 

Co and introduced earlier 


■^0,7io 


00 


E 

n—O 


1 

( 2 d)" 


E 1 

76r„(0,a:) 


® l)Apa(nxO ® 7) ■ 

n=o ' ' jer„(o,x) 


(54) 
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The coefficient -Do, 7^0 depends only on the points 0 , x but not on the path 
chosen to connect them. However, the coefficient DA,'y^o ^ priori depends on the 
path 72,0 as well, since the projected areas A^,y('jxO ® 7 ) depend on the shape of 
the closed contour 72,0 0 7 - 

3.3 Bilocal operators with a 1-hop separation 

Let us first consider the specific case where the points 0 and x are nearest 
neighbors on the lattice and the path 73,0 that connects them is the shortest 
possible, i.e. an elementary link on the cubic lattice. 

3.3.1 Coefficient Dq,* 

For the coefficient we can without loss of generality assume that the link 

72,0 is a link 72,0 = x. The random walks that connects 0 to x = 0 + x must 
have an odd length, n = 2m + 1. Simple combinatorics leads to the following 
expression, 

D - = V ^ V (2m + 1)! 

^ (2(i)2-+i ^ ni!(l + ni)!n2!2 ---nd!2 • ^ > 

By a Borel transformation similar to the one used in the section 2.6, we can 
derive the following integral representation for this coefficient® 

cOO 

Do,,= dt e-* ■ (56) 

Jo 

In 3 dimensions, this leads to a closed expression ([39]-§6.612.6), 

^ ^ 0.51638606 , (57) 

while in 4 dimensions we get 

Do,i « 0.2394671218 . (58) 


3.3.2 Coefficient D^V^ 

We must extend the discussion of the subsection 2.7 to a summation over all 
closed loops with a fixed link, of the type shown in the Hgure 5. We need to 
consider two cases: 

a. The link 72,0 connecting the two points is not in the plane )j,i> on which we 
project the areas. The archetype of this case is fiv = 12 and 72,0 = 3 (one 
elementary link in the +X3 direction). 

®In the appendix B, we derive an integral expression for this leading order coefficient 
when 0 and x are not nearest neighbors. 
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b. The link lies in the plane. For instance, = 12 and 73,0 = 1- 

The first case a is the simplest, because the fixed link does not affect in any 
way the projected area in the 12 plane. The fixed link 73,0 = 3 only alters the 
counting of the hops that are orthogonal to the 12 plane. We now have 

rji 2 ;i 2 _ _ / 2 (m + n) + l\ 

4,3 - ('2c;)2(n+m) + l I 2m + 1 ) 

X ^1 X ^ (^ 12 ( 7 ))" • (59) 

.<Ter3™+i(o.3)J ^6r3„(o.o) J 2 

Note that the number of hops in the transverse directions must be odd (we 
denote it 2m + 1). We can still use the formula (31) for the variances of the 
areas in dimension 2 (the second factor on the second line). The first factor 
on the second line counts the number of random walks of length 2 m + 1 in the 
transverse directions that connect the points 0 and 3. In d = 3 and d = 4, this 
factor is given by 



In the case b, the formula (31) must be modified in order to sum only over 
loops that start at the point 0 and end at the point 1. Let us consider random 
walks of length 2n— I, so that the length is 2n after adding a hop 1“^ to return to 
the starting point and close the loop. Consider the set of all the closed random 
walks of length 2n, over which the sum in eq. (31) is performed. After 2n — 1 
steps, these walks must be at one of the 4 four points 1, —1, 2 or —2 (since at 
the next hop they return at 0). Thus the set r2„(0,0) of these closed random 
walks can be partitioned into four subsets according to the point reached at the 
step 2n — 1. These subsets are identical up to rotations of angles multiple of 7r/2, 
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and therefore each of these subsets has the same contribution to the variance of 
the area (since the area is invariant under these rotations). Since we are now 
interested only in the random walks that reach the point 1 after 2 n — 1 steps, 
we see that eq. (31) must be replaced by^° 


E (^12(7))^ 

L76r2„-i(o.i) 


_ 1 

f2n\ 

4 

- dim 2 

V n y 


n^{n — 1 ) 
6 ( 2 n- 1 ) ■ 


Note that the combinatorial prefactor can also be written as 

2 / n._ 1 \ 2 


1 / 2 n 

din 


2 n- 1 
n 


(63) 


(64) 


which is the number of 2 -dimensional random walks of length 2 n — 1 from 0 
to 1 (see the second of eqs. (60)). From this, it is straightforward to find the 
analogue of eqs. (61) and (62) in the case where 72,0 = 1 “^ : 


^I 2 ;i 2 _ \ " (2^ + 1)! \ "_ (2n -I- 1)! _ (n -|- l)^n 

4.i d=z Q‘21+1 (^l — n)PnP{n + ly.'^ 6{2n + 1) ’ 


and 

t^ 12;12 _ (21 + 1 )! {2{l — n))l{2n + l)l (n-I-l)^n , . 

44 d=A {I — ny."^nP{n + 1)P 6{2n + 1) 

(Note that, compared to eq. (63), we have shifted n —)■ n + 1, so that the number 
of hops in the 12 plane is 2n -I- 1 instead of 2n — 1.) 

The combinatorial formulas (61), (62), (65) and ( 66 ) can be transformed into 
an integral by means of a Borel transformation. We find the following results, 
valid in d dimensions : 


and 


D 


12;12 

4,3 




4,1 


^00 

Jo 


12 d2 

1 


12d2 


(67) 

( 68 ) 


3.3.3 Behavior for long random walks 

Having in mind the infrared divergences that manifest themselves in the be¬ 
havior of the summand in for large path lengths 21, it is interesting to 

study also the summand in the newly introduced coefficients ^ and ^ . 
In the figure 6 , we plot as a function of the index I the relative differences be¬ 
tween these quantities. It appears that these differences decrease as when 

order to obtain this formula, one may also use eq. (143) from the appendix C and sum 
over 1 < ni < n, 

(2(ni-h 712 ) — 1)! nin 2 1 /2n\^ n^(n—1) 
ni!(ni — l)!n 2 !^ 3 4 \ n J 6(2n — 1) 


E 

ni+n2=n 
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d = 4 


1 



0.1 


0.01 


1 


10 


100 


Figure 6: Relative difference between the summands in the coefficients 
12 12 12 12 ^ 
D.r and D.T' as a function of the length of the random walk, in d = 4 

4,1 4,3 ® 

dimensions. 


I oo. In other words, the coefficients have the 

same behavior for large random walks. This is in fact quite intuitive, since large 
values of I correspond to long random walks, that depend very little about the 
conditions that are enforced at the endpoints of the path. 

This observation has also very practical consequences in the evaluation of 
the expectation of local operators that contain derivatives. In the particular 
example of eq. (50), it is easy to see that the second order coefficients in the 
small a expansion will always come in combinations such as or 

. The summands in these differences have a milder behavior at 
large I, and instead of having a quadratic infrared divergence, they have only a 
logarithmic divergence 



(69) 


as expected on dimensional grounds. 

3.4 Bilocal operators with a 2-hop separation 

In the discretization of (Z?i/^)*(I?^^), that appears for instance in the energy- 
momentum tensor, we also need to consider operators such as 


[x)Ul^{x)(j3c{x + fi) . 


(70) 
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When fj, = v, the two Wilson lines cancel and this operator is simply the local 
operator (j)*(j) evaluated at the point x + fl. The novel case is when ^ ^ v, for 
which the two elementary fields are separated by two hops in distinct directions. 

For the leading order coefficient JDg we can use the result derived in 
the appendix B, 

nOO 

(71) 

Jo 

For the next-to-leading order coefficient we need now to distinguish 

three cases, depending on whether the directions /i, v coincide with the directions 
1, 2 of the plane on which the areas are projected. The results for a single hop 
separation can be generalized into : 


i. fi, v ^ 1,2 : 


D 


4,i>- 


12 d 2 




(72) 


ii. ^ = 1, ^ 7^ 1, 2 : 


12;12 


1 



Ilil) > 


(73) 


iii. ly = 1, fj, = 2 : by using eq. (145) and a Borel transformation, we can 
obtain the following expression : 


D 


12;12 

4 , 1-12 
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(74) 


4 Anisotropic lattice 


4.1 Worldline representation 

Let us now consider a lattice in which the spacings are not the same in all 
directions: oi, 02 , • • • , a^. The covariant derivative squared now reads 


-D 


r 1 1 

■■ ,id) =-2 V ^ /(ii,• • • ,id) 
L aij 


r—l 


E 

r=l 


Ur{‘ • • , Zi-, •••)/(•• • 


+ 


[/-!(• •• ,V-1,---)/(••• ,ir -!,■■■) 


Let us introduce a “mean inverse squared lattice spacing”, 


(75) 


1 _ 1 1 
a2 d ^ ai ’ 

r—l ‘ 


(76) 
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and now we write 


(77) 


1 2 

n—0 

This has again the virtue of eliminating the “stationary” term in the random 
walk, since we have 




Ur(‘• • ) *r; ■■■)/(■■ ■ + 


Ur ,ir - I,--- )/(••• 


(78) 


The difference compared to the isotropic case is that the probability to make a 
hop in the direction r is modihed by the factor 


hr = 



(79) 


This factor is trivially 1 if the spacings are all equal but differs from 1 for 
an anisotropic lattice. The hopping probability in the direction r is inversely 
proportional to the squared lattice spacing in this direction (the factor hr is 
almost equal to d in the direction that has the smallest lattice spacing). Eq. (16) 
is thus generalized into 






oo 


E 

n—O 


1 

(2d)2" 


E 

7er2„(o,o) 


rt'-' 

ie-f 


tradj (w( 7 )) , 


(80) 


where 11^67 denotes the product of all the h’s collected along the closed loop. 


4.2 Geometrical interpretation 

Note that the mean squared value of the absolute distance traveled in the di¬ 
rection r is given by 

(Ml) = Ur al , (81) 

where Ur is the number of hops in the direction r. Since rir is proportional to 
hr, the product is independent of the direction r, and we have 

(Ael) = (Ael) = ■■■ = {AQ , (82) 

regardless of the values of the lattice spacings. In other words, the loops that 
enter in the formula (80) have an absolute geometrical shape which is isotropic 
(on average). On an anisotropic lattice, this isotropic distribution of shapes 
is realized by making more hops in the directions that have a smaller lattice 
spacing. 
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4.3 Leading order coefficients Cq and 


The coefficient that appears in the leading term of the continuum limit is now 
modified into 


OO .. 

E [!!'« 

n=o ^ ' 7er2„(o,o) 1^1 


(83) 


A somewhat more explicit expression can be obtained by partitioning the 2n 
hops as 2n = 2ni + - • •+2nd, where 2ni is the number of hops in the xi direction, 
2n2 the number of hops in the X 2 direction, etc... This leads to 


Co{{hr}) = 


OO 

^ (2d) 

n =0 ^ ^ 


(2n)! 

2 n 


E 


' 2nj 


niH- \-nd=n 


nl!^ • • • 


(84) 


One can easily generalize the derivation of eq. (27) in the anisotropic case. This 
leads to the following integral representation for Cq, 


Coi{hr})= dte-^ . 

Jo 


(85) 


Similarly to eq. (85), it is possible to obtain the following integral representation 
of the leading order coefficient Do,* that appears in the continuum expansion 
of composite operators involving fields separated by one lattice spacing : 


D, 


pOO _ 

,q({M)=/ 

do 


(86) 


Likewise, the anisotropic form of the formulas of the appendix B is completely 
straightforward. 

In the example of the operator of eq. (50), the leading term in the expansion 
in powers of the lattice coupling can be expressed in terms of the coefficients 
Co and Dq.*, 


m0){D^D^cl,U0)) 


f^adj (f) 
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ai 
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Do,i ~ C'o 


(87) 


Using the identities Ii = /q and {xIi{x)Y = xIo{x), and integrating by parts in 
eq. (86), we obtain 


d 

Y,h^Dor^ = d{Co-l) . 

i=l 

This identity leads to 

Do^i + 1 — Co _ 


( 88 ) 


(89) 
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( 90 ) 


When we use this property in eq. (87), it turns it into 

m0){D^DM0)) = H 

ai-'-Ud 

which is consistent with the equation of motion satisfied by the propagator : 
{K{ 0 ){D^tDf,(j))a{ 0 )) = \iin{Df,D^)f‘ {<Pl{ 0 )<Pb{x)) = -tradj (1) lim 5(a;) . 

X—>-0 rc—>-0 

(91) 

A similar identity, eq. (98), among the coefficients that appear at the next 
order ensures that this property does not receive corrections that depend on the 
background field (i.e. the dots in eq. (90) are in fact zero). 


4.4 Coefficient 

In the isotropic case, an essential ingredient in the calculation of the coefficient 
C 4 was the combinatorial formula (31), that gives the variance of the areas 
enclosed by closed random walks in 2 dimensions. Now, we need to generalize 
this formula to random walks weighted by the factors hi that depend on the 
direction of each hop. Note that in the isotropic case {hi = 1), this quantity 
has a simple generating function. 
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n—O 
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7er2„(o,o) 
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(92) 


The integral representation (85) of Cq in the anisotropic case and its comparison 
with the corresponding isotropic formula (39) suggests the following generaliza¬ 
tion of eq. (92) : 
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(4112(7))^ 
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= hih2X^ I^{2hiX) hi2h2X) . 


(93) 


It is straightforward to see that this formula is equivalent to the identity (142) 
listed in the appendix C. In a sense, eq. (142) can be viewed as a “fine grained” 
version of eq. (31), that retains the information about the number of hops in 
each direction. We present a proof of this identity in a separate paper [38]. 

With the help of eq. (93), it is immediate to obtain the integral representation 
of the coefficient in the anisotropic case. 


^12;12 


{{hr}) 




dt 


[ n ^o(¥)] Ii{h^)h{hf) 


*#1,2 


(94) 


This coefficient will enter in the short distance expansion via the combination 


J,1U2 


2 (0)A,12(0) ^ a4F,i2(0)Ff (0) , 


(95) 
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where the unwritten factors are dimensionless numbers. Therefore, the terms in 
are accompanied by the fourth power of the “average lattice spacing” a, with 
a numerical prefactor that depends on the ratios aijaj of the lattice spacings in 
the various directions. Thus, if we take the short distance limit while keeping 
these ratios hxed, only the factor decreases, while the numerical prefactor 
stays constant. 


4.5 Coefficients and 

Likewise, eqs. (67) and (68) can be generalized into 


nl2;12/r, IN hih2 


dt e * 


n ^o(^)]/i(^)/i(^)/i(^), 

i#l,2.3 

(96) 


and 


= ^ j^dte-U^ [n^o(^)] /i(^). (97) 


Using these formulas, we can also prove the identity 
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Y.h,D 


12;12 

4,i 


i=l 


= dCi 


(98) 


that generalizes eq. (89) to the coefficients that appear in the next order of the 
expansion. 

It is also easy to generalize to the case an anisotropic lattice the coefficients 
with a separation of two hops, given in eqs.(72), (73) and (74) in the 
isotropic case. One obtains : 

i. 7^ 1, 2 : 
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(99) 


ii. ^ = 1, ^ 7^ 1, 2 : 
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dt ( 
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n ^o( 

2^1,2,/i 


hi t 


), 

( 100 ) 


word of caution is necessary here: the l.h.s. and r.h.s. of this equation contain infrared 
divergences that should be properly regularized, for instance by the introduction of a mass 
(see the section 2.8). 
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iii. V = 1, ^ = 2 : 


pjl2;12 

■^ 4 , 1-12 


dt 


hih2t'^ 

L 2d? 


^o(^)/o(T) 


n ^o(^). ( 101 ) 


1 , 41,2 


4.6 Limit of extreme anisotropy 

An extreme case of anisotropy is to have one lattice spacing much smaller than 
the others. This could correspond to the case where one of the coordinates (e.g. 
the time) is treated as a continuous variable while the other are discretized with 
a small but finite lattice spacing. Let us therefore assume that 


Cd < ai, • • • , Od-i0 . (102) 

In order to get a sense of the behavior of the coefhcients in this limit, let 
us consider the coefficient Co, whose integral representation (85) can easily be 
studied numerically. In this numerical study, we consider the dimension d = 4, 
and the remaining 3 lattice spacings are assumed to be all equal. As shown in 
the figure 7, we see that Cq diverges as in this limit. In this limit, we have 



84 / ai 


Figure 7: Behavior of the leading coefficient Cq as a function of the ratio of 
lattice spacings 04 / 01 , in 4 dimensions. 
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404 


a? ^ 4 o 4 

O1O2O3O4 O1O2O3 


(103) 
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so that the leading term of eq. (80) in fact becomes independent of the smallest 
lattice spacing. This result is consistent with the result of the appendix F, where 
we use from the start a continuous time variable. 

It is in fact possible to understand this limit analytically, starting from the 
integral representation of the coefficients that appear in the small lattice spacing 
expansion. Generically, these coefficients involve integrals of the form, 


An-,{Si} 


dt i 




(104) 


where the exponent n is 0 or 2 and the indices Si are 0 or 1, in all the examples 
we have encountered so far. In this limit, we have 


hd = d 



{i < d) . 


(105) 


Therefore, the argument of the last Bessel function is much larger than the 
arguments of the first d — 1 Bessel functions. Let us first define a rescaled 
integration variable by 


-(- 0 ) 

, dsal 
di + dsal ’ 

(106) 

ds=d-l , 

1 1 1 

a? " rfs ^ df ■ 

(107) 


where we denote 


Using hit/d = SL^r/^dsaf) and the asymptotic expansion of the modified Bessel 
function. 


Isiz) = 


1 


z->-+oo 


+ C>(z-3/2) 


(108) 


we arrive at 


-Arj.- 


1 / a: 


(sij) I (100) 


n;{Si} 


o-O 


where we define his = A^lo-i- Note that this limiting value does not depend on 
the index Sd of the Bessel function associated to the temporal direction. 

When n = 0 (i.e. for the coefficients Co and Dq), this integral behaves as 
which is cancelled by the behavior of the overall prefactor sd/{ai ■ ■ - ad)- 
Therefore, the leading order terms have a finite limit^^ when ad —> 0. In the 


we take ai = ■ ■ • = a^_i ^ a^, then the leading term of (0)</)a(0)) reads 


{ra{0)M0)) 




tl-adj (1) 1 r°° 

2af“^ '/dmTs Jo Vr 



0.4553440518--- if d^-Z 


in perfect agreement with the figure 7. 
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next-to-leading order coefficients and the integral (109) appear 

with the exponent n = 2 and therefore it behaves as ■ One of the powers of 
is cancelled by the prefactor a ?/(ai • • • ad)- Moreover, the coefficients 
and also contain a prefactor hihj, 


hihj 


9 9 "> 

ad^o afa^ 

*' j 


( 110 ) 


thereby canceling the remaining factor from eq. (109). Therefore, the next- 
to-leading order coefficients are also finite in the limit —)■ 0 (aside from the 
possible infrared divergences discussed previously). 


4.7 Energy-momentum tensor 

Thanks to the results of the previous sections, we can write the short distance 
expansion of the energy momentum tensor. In this section, we present expres¬ 
sions for its diagonal components (its off-diagonal components can be treated 
similarly, but have a slightly more complicated structure). These diagonal com¬ 
ponents read : 

T" = (A<^)*(A0) - . (Ill) 

3¥=i 

Therefore, we need the expression of the expectation value of the operator 
{Di(f))* (Diip) (not summed over i). The lattice version of this operator is 




r (o)<^(o)+r ml {o)m - r mi^) 

( 112 ) 


and the worldline representation of its 1-loop expectation value reads 


2 1 1 

i n—0 ^ ' 


[n^*] (w( 7 )) 

7er„(o.o) ie-y 

E [n foj tradj (W(7)) 

7Gr„{f,z) ^^7 

■ E [n^ tradj (t//(0)W(7)) 

7 er„(i.o) ie-y 


E 

76r„(0.i) t&'y 


t^adj (t/70)W(7)) . (113) 
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Using results derived earlier in this paper, we can write the following short 
distance expansion for this quantity : 




2 d Hr®*' 


tradj ( 1 ) {Co — -Do,i) 


-T E +• • • • (114) 


fJ-Ki' 


Using eq. (Ill), as well as eqs. (89) and (98), we obtain the following expansion 
for the diagonal components of the energy-momentum tensor. 




dUr 


2 a: 


I tradj ( 1 ) ^^2 (^0 -Do,i) ^ 2 ) 
2 ^ alaf{F^\0)y (Cf 


1 k<l 
k 


9 \ ^ 2 2 { T7>ik (c\\\‘^ {/-iik\ik j^ik;ik\ 

W) {C4 -D^: ) 


k^i 


(115) 


where we have explicitly separated the terms of order depending on whether 
one of the indices carried by the field strength is i or not. Integral expressions 
in terms of modified Bessel functions for all the coefficients that appear in this 
formula, for arbitrary lattice spacings, can be found in the previous sections of 
this paper. By using eqs. (89) and (98), we can see that 


= ( 2 -d)*^^l^-f ••• , (116) 

7^1 

where the dots are UV finite terms in four dimensions. 

Eq. (115) provides the final answer for generic lattice spacings. If in addition 
we assume that the temporal spacing is much smaller than the others, we 
can use 




ad-^o aj 


+ 


f 2 a 2 


\Trdsal 


1/2 


dr 




2 ( , _ 1 
n^. ( ad->0 a^- ’ 


(117) 


d 


respectively ior i ^ d and i = d. For the next-to-leading order coefficients, 
we can also use eq. (109) to obtain their limiting value if z 7 ^ d. If z = d, it 
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is necessary to go one order further in the asymptotic expansion of the Bessel 
function of argument hdt/d. In this case, we can use the following limit 



The terms in eq. (117) lead to contributions proportional to (ai • • • ad)~^- 
It is easy to check that all the other contributions to T" become independent 
of ad in the limit —>■ 0. 


5 Conclusions 

In this paper, we have applied a discrete version of the worldline formalism in 
order to obtain expressions for 1-loop expectation values in a lattice scalar field 
theory, in the presence of a non-Abelian gauge background. In this framework, 
2 -point correlators are expressed as sums over all the random walks that connect 
their endpoints (or closed loops in the specific case of local operators). This rep¬ 
resentation renders the ultraviolet and infrared behaviors of these expectation 
values very intuitive: ultraviolet divergences are encoded in the contribution of 
very short random walks, while the infrared behavior arises from the statistics 
of long random walks that explore large regions of space-time. 

Moreover, it is straightforward to take the limit of small lattice spacing in 
these worldline expressions. This gives an expansion in powers of the background 
field strength, whose coefficients are sums of all the closed random walks on the 
lattice, thereby reducing the calculation of the coefficients to a combinatorial 
problem. The leading coefficient is merely counting these random walks, and 
is therefore very easy to obtain. The next-to-leading order coefficient is related 
to the variance of the areas (projected on a plane) enclosed by these random 
walks. The relevant combinatorial formulas can be found in ref. [37]. 

When one considers a lattice with anisotropic lattice spacings, these random 
walks are further weighted by factors that count the number of hops in each 
direction. In order to obtain simple expressions for the coefficients of the ex¬ 
pansion in this case, we had to conjecture some generalizations of the formulas 
of ref. [37]. The proof of these formulas is given in a separate paper [38]. 

Using a Borel transformation, all the coefficients that appear in the short dis¬ 
tance expansion of these correlators can be rewritten as I-dimensional integrals. 
These can be easily evaluated numerically, and they are also quite convenient 
in order to study analytically the limit where one lattice spacing becomes much 
smaller than the others. 

This work can be extended in several directions. One of them is to depart 
from a scalar field theory, in order to study quantities where a spin 1/2 fermion 
or a spin I gauge boson circulates in the loop. In the continuum theory, these 
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extensions are well known. Another -considerably more difficult- extension 
would be to use the worldline formalism in order to go beyond one loop. Also, 
at the moment, it is unclear whether one could modify this formalism to handle 
a Minkowskian time (see the discussion in the appendix F). 
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A Closed random walks on a periodic lattice 

In the main body of this paper, all the combinatorial formulas that we have 
derived count random walks on an infinite cubic lattice. Since the lattice is 
infinite, there is no need to specify boundary conditions. This is the appropri¬ 
ate setting when one considers the limit a —)■ 0 at fixed physical volume: the 
constant volume is achieved by letting the number of lattice points go to infinity 
(as a“^) in each direction. In this appendix, we consider a different limit, where 
we keep fixed the number of lattice points as we take the limit a —>■ 0. We now 
assume periodic boundary conditions. 

We have introduced in eq. (9) a sequence of functions defined on the 
lattice. If the first of these functions is 


Po{x) = 4,0 , 


(119) 


then Pn (x) is the probability that a random walk starting at the point 0 reaches 
the point x after n steps^^. With this choice of initial condition, P„(0) is the 
probability that the random walk makes a closed loop of base point 0 in n steps. 

The first step is to find the spectrum of the linear operator that maps 
to Pn+i- By using eq. (13), one finds that the eigenfunctions and eigenvalues of 
this operator are the plane waves 


'/’fc.d(^) = , {it ■ x}d = kixi - \-kdXd, (120) 

and that the associated eigenvalue is 



^^Using eq. (9) and “integrating by parts”, it is easy to check that the quantity 

yelattice 


is conserved as n increases. 
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where N is the number of lattice spacings in each direction (for simplicity, we 
assume the same number of spacings in all the directions). The initial condition 
of the iteration can be written as 

{k}d 


where {k}d denotes all the d-uplets of integers in the range [0,N — 1]. After n 
iterations this distribution has become 

{k}d 

In order to obtain the probability for a random loop of length n to be closed, 
we simply evaluate this at a; = 0, 

E ■ (124) 

{k}d 


The sum over n of the probability of closed random walks of length n therefore 
reads 

Ei'40) = Ty;^, (125) 

n>0 k,d 

where we recognize the lattice expression for a 1-loop tadpole. 

Let us now consider in more detail eq. (124). We can make the right hand 
side more explicit by writing 


1 


k,d 

{k}d 


Nd{2dy 




" ' e ^ 


OjirLi _27vLi^" 

-1_ 5^*44 „ + e " 


{fc}d 


1 V ^ V ^ Tl\ 

N<^{2d)'^ ^ ^ nilpil ■ ■ ■ Udlpd'- 

{k}d "i+i’i+ ■ 


+ "d+P£i = 


X e « • • • e « 


E 


(2(i)” ni!pi!---nd!pd! 

+"d+Ed=" 


I I (^ni-pi,0[Af] ■ ■ ■ (^nd-pd.0[A(] ) 


(126) 


where the symbol (5p,o[7V] 
N : 


is a Kronecker symbol that defines an equality modulo 


0p.0[A(] 


- E'^p- 


iN ■ 


(127) 


2 G S 

The essential difference between an infinite lattice and a finite lattice with pe¬ 
riodic boundary conditions comes from here. 
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If we take the limit iV —>■ oo, then we have 


Op,0[Af] 


= 

N—>oo 


'p,0 1 


(128) 


and eq. (126) is equivalent to 


1 

w 


{k}d 


1 

(2d)2™ 

0 


E 

niH- \-nd=m 


(2m)! 

ni!^ • • • 


(n 


2m even) 

, (129) 

(n odd) 


which is nothing but the combinatorial expression for the coefficient Cq on an 
infinite cubic lattice. 

In order to illustrate explicitly the difference on a Hnite periodic lattice, let 
us consider a 1-dimensional lattice of N sites. Eq. (126) becomes 


ivE^I.i- 2"Ep!(^_p)!'^2p.„[Ar] • 
{fe}i p=o^^ 


(130) 


For N even, this forces the random walk to have an even length n = 2m, and 
we get 

— = 

TV ^ N even 

{fell 


where the sum in the right hand side stops when the argument of the second 
factorial in the denominator becomes negative. The first term of this formula is 
identical to the result on an infinite lattice. The additional terms corresponds 
to paths that wrap around the periodic lattice, with winding number ±1 for 
the second term, ±2 for the third term, etc... These terms can be viewed as 
finite size corrections, since they explicitly depend on N. When the lattice size 
N is odd, the situation is even more complicated because random walks of odd 
length are also permitted. For n = 2m, we get 


22” 


(2m)! 


-k2- 


(2m)! 


L m!2 (^ui + y)K''tT'~ 


+2 


(2m)! 


(m - 


^)!(m-^)! 


-k • ■ 


, (131) 


— V 

JY k,l jv odd 

{fell 
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- (2m)! ^ ^ (2m)! 

- m!2 (m + N)!(m 


N)! 


o (2m)! 

(m + 2N)!(m-2N)! ’ 


(132) 


and for odd lengths n = 2m + 1 we have 


Af ^ 


2m+l 


1 


N 


(2m -k 1)! 


{fell 


fc.i AT odd 22-L (to+ i^)!(m+i^)! 

■ r, (2m-kl)! 

(m+i^)!(m+i^)! 


.(133) 
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The situation in d dimensions follows the same pattern, but explicit formulas 
become quite cumbersome. In particular, since the random walk takes place on 
a d-dimensional torus instead of a ring, the “winding number” is now a d-uplet 
(wi, • • • jWd) G For a given length n of random walk, the term of winding 
number (0, • • • ,0) is identical to the result on an infinite lattice, but in addition 
one must sum over all the possible windings allowed for a given length n of the 
random walk. 


B Coefficient for arbitrary separations 


In this appendix, we generalize the eqs. (55)-(58) to the case where the points 
0 and X are separated by more than one hop. For dehniteness, let us assume 
that the point x is 

X = xii + • • • + Xdd . (134) 

The minimal length of a path connecting 0 to a; is: 

A = xi + ■ ■ ■ + Xd . (135) 


Firstly, let us notice that the parity of the length of the paths jxo that connect 
the two points depends solely on the endpoints 0 and x. This parity is that of 
the number A. One can first obtain a combinatorial expression for Dg -y^^, 


D, 


0,7x0 


= E 


E 


ni 


+ Xl)'. ■ ■ ■ Ud'-iUd + X,)! 


(136) 


By a Borel transformation, this can be turned into 


D 


0,7x0 


COO 

/ dt e"‘ Bd(^) , 
-'0 


with 


Bd{x) = 


2=1 


E 

.p—O 


^2p-\-Xi 


pl(p + Xi)l 


= Y[Ix,{2x) . 




This leads to the following integral representation of 


D, 


0,7a: 


dte ‘ 47 


2 = 1 


(137) 

(138) 


(139) 


In 3 dimensions, this leads to a closed expression (see ref. [39]-§6.612.6) 


with 


= 3 ( rig + 

^2,1 
-1" ’"3 1 ) 

(140) 

d—3 y 

7r2g J 

>/3 ~ 1 p2 1 
967r3 ''24 


(141) 


and where the coefficients ri^ 2,3 are listed in the table 2 as a function of a;i, 2,3 for 
small separations. In 4 dimensions, the integral of eq. (139) must be evaluated 
numerically. Some values for small separations are listed in the table 3. 
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a;ia;2X3 

ri 

r2 

rs 

000 

1 

0 

0 

100 

1 

0 

-1/3 

no 

5/12 

-1/2 

0 

200 

10/3 

-2 

-2 

111 

-1/8 

3/4 

0 

210 

3/8 

-9/4 

1/3 

300 

35/2 

21 

-13 


Table 2: Values of the coefficients ri^ 2,3 for separations up to three hops (from 
[39]-§6.612.6). These values are invariant under permutations of a:i^ 2 , 3 - For 
larger separations, see [39]-§6.612.6. 


Xia;2a:3X4 

-^0,a;i,2,3,4 

0000 

1.239467122 

1000 

0.2394671218 

1100 

0.1017176302 

2000 

0.06596407193 

1110 

0.06187238110 

2100 

0.04365863661 

3000 

0.02629363394 


Table 3 : Values of the leading order coefficient -Do,7^0 fo ^ dimensions, for 
separations up to three hops. These values are invariant under permutations of 

2^1, 2 , 3 , 4 - 
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C Moments of the distribution of areas in d = 2 

C.l Variance of the areas of closed loops with fixed sec¬ 
tions 

In the case of bilocal operators, the next-to-leading order coefficients 
require that we evaluate the variance of the areas enclosed by the paths 
where 7 is a 2-dimensional random walk from the point 0 to the point x. In 
other words, 7 ® is a closed path, but only the section 7 is random, while 
the section 73,0 is held fixed. 



Figure 8 : Left: closed random walk from 0 to 0. Right: closed random walk 
with a fixed section of length 1 between 0 and x (shown in green). 




Figure 9: Closed random walks with a fixed section of length 2 (shown in green). 
Left: a; = 2 i. Right: a: = i -I- 2. 

Having in mind the case of an anisotropic lattice, it is also useful to have 
formulas that keep track of the number of hops in the xi and X 2 directions sep¬ 
arately. All the formulas listed in this appendix have been checked numerically 
by an exhaustive sum over all random walks of length 2n < 20. In ref. [38], 
we present a proof of all the formulas of the next subsection, where the two 
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endpoints are identical. The same method can be used to study the case where 
the endpoints are distinct, although we have not yet done so. 

In the following, we give formulas for the four cases illustrated in the figures 
8 and 9, which covers all the situations where 0 and x are separated by two hops 
at most. The simplest situation is when the point a; = 0 and the path is the 
null path (figure 8 , left). In this case, we have^^ 


(Area ( 7 ))^ 

7er„j,„2(o.o) 


(2(ni + 712 ))! nin 2 
ni!^n2!^ 3 


(142) 


(We recall that (0, x) is the set of all the paths that connect the point 0 to 

the point x and have ni hops in the direction +xi and n 2 hops in the direction 
+X2.) 

Next, let us consider the case where x = xi and (i-®- ^ single hop 

in the —xi direction - see the figure 8 , right). This case leads to 


(Area( 7 (g)i ^))^ 

(Oil) 


(2(ni + 712 ) - 1)! 771772 
77i!(77i — 1)1772!^ 3 


(143) 


For separations of two hops, there are two possibilities, shown in the figure 9. 
When the point x is located at a: = 2 1 and the path connecting to 0 is made of 
two horizontal hops (figure 9, left), we have 


(Area( 7 (g)i ^))^ 

7 er„j ,„2 (o, 2 i) 


(2(771 + 772 - 1))! (771 + 1)772 
77i! (771 — 2)1772!^ 3 


(144) 


When the point a; is at x = 1 + ^ and is connected to the point 0 by the path 
7 a ;0 = 2 “^ 1 “^, we obtain 


y^ (Area ( 7 ® 2 ^))^ 

( 0 , 1 + 2 ) 


(2(771 + 772 - 1))! / 771772 1 \ 

77i!(77i — 1)!772!(772 — 1)! V 3 6 / 

(145) 


C.2 Higher moments (for a; = 0 and 7^,0 = 1) 


In the case of closed paths 7 , it is also possible to obtain by empirical observation 
a simple formula for the moment of order 4 of the area, 


(Area ( 7 ))^ = 

7er„i.„7o,o) 


4 (2(771 + 772))! 77 1 772 


77i!2772!^ 


15 


(7771772 - (771 + 772 )) , (146) 


can check that 


(2n)! ni(n —ni) (2n)!^ n^(n—1) 

— ni)\‘^ 3 6(2n — 1) ’ 

which is nothing but the isotropic result for the variance of the areas in the set of all the 
closed random walks of lengths 2n (see [37], eqs. (1.4)—(1.5)). 
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that leads to the formula (1.6) of ref. [37] after summing over 0 < ni < n 
(with n 2 = n — ni). Interestingly, the formulas for the moments are somewhat 
simpler if one considers only random paths with fixed numbers of hops in the 
-\-Xi and +X 2 directions, rather that all the random paths that have a fixed 
length. Indeed, besides the expected combinatorial factor, these formulas seem 
to involve a polynomial in ni^ 2 , instead of a rational fraction in n = ni + n 2 . 

Similarly, we have obtained the following expressions for the moments of 
order 6 , 8 and 10 : 

/A / n \6 ( 2 (ni+n 2 ))! UiUa . , >2 

— I5nin2(ni + 712) + 2(ni + 712)^ — (ni + 712)) , (147) 

7 a / a\8 ( 2 (ni+ 772))! 771772/ , .3 

-134(771772)^(771 + 772) + 5377i772(77i + 772)^ - 6(77i + 77^) 

-4077i772(77i + 772) + 8(77i + 772)^ - 3(77i + 772)) , (148) 


E 

7Gr.n,^ ^^2 ( 0 ) 0 ) 


E 

(0,0) 


7a / a \10 ( 2(771 +772))! 771772 / , ,4 

(Area ( 7 )) = -^ (2555 (771772) 

—4778(771772)^(771 + 772) + 3745(771772)^(771 + 772)^ 

— 5290(771772)^(771 + 772) — 1282771772 (77i + 77^) 

+ 1918771772 (771 + 772)^ + 120(771 ~ ^^2)^ 

— 1403771772 (771 + 772) — 300(771 + ^2) 

+ 270(771 + 772)2 - 85(771 + 772 )) . ( 149 ) 

These observations suggest that we have in general 

E (Area (7))"'' = V2k{ni,n2) , ( 150 ) 

7 er„^.„,(o.o) 2 - 

where 7 ^ 2 fc( 77 i, 772) is a polynomial of degree 2 k, symmetric in (771,772), propor¬ 
tional to 77 1 772 ,with rational coefficients^®. The leading term of this polynomial 
is of the form 0^(771772)^”^, and the first five terms obtained above are consistent 
with (defining cq = 1 ) 

<-) 

^^Such a polynomial has k?‘ independent coefficients. If we evaluate the areas of all the 
closed paths with ni,n 2 such that ni + n 2 < n, we get independent constraints n = 2p 
and p(p+ 1) independent constraints if n = 2p+ 1. Therefore, in order to determine uniquely 
the polynomial 'p 2 ki it is sufficient to consider all the closed paths up to the length 2n = 4/c. 
This is why one can obtain all the moments up to 2k = 10 by considering all the paths up to 
the length 20. 


E 

76r„j.„2(0,0) 
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in agreement^® with the eq. (1.7) of ref. [37]. In a separate paper, ref. [38], we 
prove eq. (150) that gives the general structure of the moments. The approach 
used in this paper also provides an algorithm for calculating explicit the poly¬ 
nomial V 2 k for small values of k, and we have checked the formulas (142), (146), 
(147), (148) and (149). 

D Link with the almost-Mathieu operator 

One can relate the statistics of closed random walks on H? to Euclidean lattice 
scalar QED on a 2-dimensional lattice by means of the discrete worldline for¬ 
malism. In order to keep track of the areas enclosed by these random walks, 
one should add a magnetic field transverse to the plane in which the charged 
scalar particles live. Moreover, since we would like to distinguish the random 
walks according to the number of hops they make in the x and y directions, we 
need a rectangular lattice, with distinct lattice spacings Oi and 02 in the two 
directions. In order to write explicit expressions, we must choose a gauge for 
the magnetic field. A convenient choice is 

^i(*.j) = 0 , A2{i,j) = Bx = Baii, (152) 

where {i,j) are the integers that label the position of a point on the lattice. In 
the lattice formulation, this background electromagnetic field is implemented 
by gauge links^^ 


C/i(i,j) = l , [/2(i,j)=e*^“^“iV (153) 


In the following, we denote 2'k/ 3 = Baia 2 the magnetic flux through an elemen¬ 
tary plaquette of the lattice. With this gauge choice, the action of the inverse 
propagator on a test function reads 


(G-V)*o 





j+i + 


j-i 


j' 


a 


2 

2 


(154) 

Moreover, the worldline representation of the propagator at equal points (the 
choice of the point 0 is arbitrary, and irrelevant since the magnetic field is 
homogeneous) reads 


a 




4 X! 42(ni-Hn2) 

ni,n2=0 7er„j,„2(0,0) i&l 


H [11''' 


^^For simplicity, the electrical charge is taken to be e = 1. 


(155) 


order to establish this connection, it is sufficient to notice that 

(2n)! fe (2n)!2 l fn^k , ^ . ] 

> ——-^ n, = -— — -e subleading terms m n . 

ni!2('n-nil!2 ^ n!4 LV2/ ^ 1 
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where we have defined 


2 



hi,2 = 


‘‘ 1,2 


(156) 


(Note that hi + h 2 = 2.) r„j_„ 2 ( 0 , 0 ) is the set of the closed paths (from (0,0) 
to ( 0 , 0 )) drawn on the lattice, that have exactly ni hops in the +x direction 
and 712 hops in the +y direction (and therefore also ni ,2 hops in the —x and 
—y directions). The product of the link variables encountered along the closed 
path is also the exponential of the magnetic flux, 


J|[/^ = (157) 


where Area ( 7 ) is the algebraic area enclosed by the path 7 . Therefore, we have 


G(0,0) 



711 , 712—0 


42 (ni+n 2 ) 


E 


i 27 r /3 Area (7) 


,, 2 ( 0 ^ 0 ) 


(158) 


and one can view the diagonal elements of the propagator in a magnetic field 
as a generating function for the distribution of the areas of closed loops on the 
lattice. 

The links exp(z27r/3i) depend only on the i coordinate. Therefore, one can 
perform a Fourier transform on the j coordinate. Let us define : 

Af-l 

(t>i,k = X! , (159) 

1=0 


where N is the number of lattice spacings in the j direction. The conjugate 
index k is also an integer in the range [0, IV — 1]. The inverse Fourier transform 
reads 

1 

= ]y E 

k=0 

If we consider an infinitely large lattice, N —>■ +oo, one can use a continuous 
momentum variable ly = 27rfc/iV, so that the above equation becomes 

^ A,, ■ (161) 

By inserting this equation into eq. (154), we obtain 






(162) 


where H^h) is known as the (anisotropic) almost-Mathieu operator : 

= hi{5i,ii+i+5i,i'-i)+2h2Cos{2TTj3i + iy) 5i,i' 2 'k5{v—v') . (163) 
V / L J 
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This Hamiltonian has been the subject of many studies, both for its practical 
interest in models of the quantum Hall effect [40, 41, 42, 43], and for its intrinsic 
mathematical interest [44, 45, 46, 47, 48, 49, 50, 51] as an example of quasi- 
periodic Hamiltonian (when j3 is an irrational number). 

From eq. (162), we get that 

2 ^ 

n—0 

Since the external magnetic field is homogeneous, the diagonal elements of the 
propagator can be obtained by taking the trace (divided by the number of lattice 
sites in the i direction). Thus, by comparing eqs. (158) and (164) we recover 
the well-known connection[52, 53, 54] between the trace of and the 

distribution of the areas of the closed random walks of length n. 





gi2,r/3Area(7) ^ lim ^ ^ 

N^oo JM Jq 27r L 


E {^0*(O)0a(O)) from lattice perturbation theory 

In this appendix, we calculate ((()*(a;)^a(a^)) in lattice perturbation theory in 
order to make the connection with our results obtained in the lattice worldline 
representation. In lattice perturbation theory (see ref. [18] for a comprehensive 
review), the expectation value of any operator O can be calculated perturba- 
tively according to 

In our case, the action reads 

S[(j)\(j>] = So[(j>\(j)] + , (167) 

with 



So[4’\(p] = -Va'^ (1)1'^ ^ [<p{x + Urf) + (j){x - arf) - 2(j){x)] ,(168) 

X r—1 ‘ 

1 p 

5'i„t[^^(('] = -Va (f^ixYUrix) - l)(j){x + arf) 

1 ■ 

X r—1 

+(j)'^ {x + arf){U~^{x)—I)<j>{x) , (169) 
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and Va = Oi- In the limit of an infinitely large lattice, the free propagator 

r—0 

reads 


Goix,y) = 


d ^JL. 


n 


dpr 


,ipr{Xr--yr 


1 


2 E ^[1 - cos(pia/)] 
1=1 ‘ 


Va 


n 


dp, 

27r 




1 


2 E ^(1 - cospi) 
1=1 ‘ 


(170) 


G({pd) 

with Xr = mrar,yr = rirar- Since it is a Green’s function of the discrete 
Euclidean d’Alembertian operator, it satisfies the identity 

^' 1 1 
- X! ^ [^o(a; + arf, y) + Gq{x - a^r, y) - 2 Gq{x, y)] = —. (171) 


1 

r-=l ' 


Let us now calculate the first two terms of the expansion of ((j)’^{0)(j)a{0)} in 
powers of the external field. 




(<).:(0)<^a(0)) = A y + A y + 


t ; 


' ■ = Iq + II + I 2 


The term of order zero is given by : 


(172) 


/o = A y = 


:V_tGdj(l) r d‘^P 


\ y 2dVa 1 


E O/ cos Pi 


tradj(l) 

2dVa 

a^ tradj(l) 




2dV„ 


dt e 


dt i 


1=1 

d’^p ( hrt 


a Jo 


(27r)'’* 

nw¥) 


I ^ ^ /Lro 

exp 2^ — cos Pi 


a=i 


(173) 


r—1 


which is equivalent to eqs. (80) and (85). 

Next, we calculate the corrections to (<))* (0 )())q( 0)) of order 0{g^a^). By 
using the power series expansion of the link variable, 

„2 

Ur{x) = g-’-ja-rArix) _ j _ iga^Ar{x) — —alA^{x) + • • • , (174) 
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and keeping only terms of 0 {g^), we find 

^2 r jdi d 


h+l 2 = ^ 

with 


y / E A:{kM:ik 2 )Irs{h,k 2 ) + 0{g^) , (175) 


Irs{ki, ^ 2 ) = [ - 77 ^ G{{pi})\^—d{{pi + kiiai})G{{pi - fezo/}) 

Va (2Tr)d lartts 

X _giPr+iPs-ik2Bas _^g-i{Pr+kirar)+iiPs-k2Ba.s) _^-iiPr+kir.ar.)-ips'^ 

-^G({pz + (Aziz + fc2z)az})(e*P’' + . (176) 

The terms of order g^a^ can be obtained by performing a series expansion of 
Irs in terms of fci,2a. It is easy to show that 


with 


Irsikljk 2 ) — klrs T l 2 ,lrs T ^ 2 , 2 rs T l2,3rs + l 2 , 4 rs + G{a^) , (177) 


J- _ Srs dp ^ J l2 2 G{{pi}) 

klrs — 2y -E J (271)^^ {kli + k 2 i) Oj G({pz}) ^^2 COSPr 


/2,lrs - y E 


2 aiaj 


Va ( 27 r)‘^ Qrtts 

1,3 = 1 


{kiikij + k 2 ik 2 j) G^({p;}) 


d^G{{pi}) . 

^ O O -sm Pr sm ps , 

UPiUPj 


I2,2ts - - y E 


d'^p OiO, 


K 2^ J_^ { 2 ttY Gra, 

k-,J - J- 


kiik2jG{{pi}) 


dG{{pi}) dG{{pi}) . 

X----- smpr smps , 

upi ^Pj 


l2,3rs — — 


2 4d 


^2.r r.,. dG{{pi}) 


Va J-. { 2 ^Y 


G{{pi}) (kls-k2s)klr 


dps 


sm Ps COS Pr 


n , \i 9G{{pi}) 

-(kir - k 2 r)k 2 s -X-cosp^sinp^ 

upr 


-^2,4rs — 


\ Yid 


G^({pz}) |^drs(fci^ + A:|r)sin^Pr 

+kirk2s (cOSPr COSPs + sin^p^d^r) 1 -(178) 


K 7_. (27rY 
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In order to simplify the above integrals, we need the following identities 


dG{{pi}) 

dpr 

d^G{{pi}) 

dprdps 


-^sinp,. G'^{{pi}) 


sinpr sinps G^{{pi}) -^ cospr G'^{{pi})6r 


= ‘2G-\{pi}) 


dGiipi}) dG{{pi}) 2 


dpr 


dps 


-^ COSPr G\{pi})5r 


(179) 


Moreover, the terms that are total derivatives with respect to any component 
of p in the above integrands always give vanishing contributions. After some 
algebra we find 

hrs = ^ [ 7^^ [ '^{ku + k2if COSp, COSPr 

+ 2{kir + k2r)‘^ Pr] , 


1 df'p r 

^2,lrs ~ ry-tz I 7r> w7 G ({t*/}) ‘2{kirkis + k2rk2s) ^*^^Pr ^^^Ps 

314 J-TT L 

d 

— {k\j. + k2r)5rs Sin^Pr — ^(^li + ^2i) COSPi COS Pr Sr 

i=l 


l2.2rs — — 




G\{pi}) {kirk 2s + kisk2r) COSPr COSPs 


3K (27r)rf 

d 

— 3kirk2rSrs Sin^ Pr + fclifc2i COSpi COSPrSr 


2 d’^p '' r 

I2.3rs=-Tprr ley w G^{{pi}) {kiskir + kiskir-2kirk2s)cOSPsCOSPr 
3l47-7r(27I')‘‘ L 

— {klr + klr — 2kirk2r)Srs S\V? Pr ■ (180) 

As a result, up to 0{g^a^), we have 

1 d'^p r 4^' 

Irs{klyk 2 ) = Tprr i ,ry Srsy^kiik 2 iCOSPiCOSPr 

314 J-y, (27r)“ L ^ 

— kisk2r COSPr COSPs] , (181) 
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and 


i-TT Jd 


h+h = - 


^^Yl[9rA'^(x) - dsA'^{x)]^J G^i{pi}) COSPrCOSPs 


192d3Va 


^[a^A“(a;) - dsA’i.{x)f 


r^s 


pOO _ 

X / dtfe-*h{^)h{!^) n ^o(^), 

-'0 , . 


i^r,s 


where we have used the integral 

1 _ 1 
a;3 2 


dt t e 


2 ^ — tx 


(182) 


In summary, we have 


(C(O)^a(O)) = 


' tradj(l) 


2dV„ 


a JO 


dte-^ 


r—1 


g^a6 

192^314 


Y.idrA’iix) - 


r^s 


pOO _ 

X / dtfe-*hi^)h{^) n ^0(^) + O(gV). 
JO „ 

(183) 


i^r,s 


This result is identical to the formulas obtained in the worldline formalism for 
this expectation value (see the eqs. (80) and (94)). Therefore, one can view this 
alternate derivation as an indirect analytical proof of the formula (142), that we 
have used in obtaining the next-to-leading order coefficient in the anisotropic 
case. In principle, one could apply the same technique in order to prove all the 
other formulas conjectured in the appendix C, although this appears to be a 
rather intricate task. 


F Continuous time and discrete space 

The extreme case of anisotropy is realized when one of the coordinates (time 
in our discussion) is treated as a continuous variable, while the others remain 
discretized on a lattice of spacing a. In order to simplify the treatment of the 
background field, it is useful to assume that the temporal gauge is used, Aq = 0. 
Therefore, the inverse propagator is 

dg 

D^=9ood^-T.^t^^ , (184) 

2=1 
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where are the forward and backward discrete covariant derivatives on a 
grid of lattice spacing a, and ds = d — 1 the number of spatial dimensions. 
At this point, we have also kept undetermined the 00 component of the metric 
tensor, in order to discuss later the difference in this formalism between the 
Minkowski and the Euclidean metric. 

For the purpose of this discussion, we can first ignore the background field 
completely, and reintroduce it later via the Wilson loop made of the gauge links 
accumulated along the random walk. Let us recall that 


V+V-/(x) = 


f{x + i)+fix-i) 2/(x) 


It is convenient to write the inverse propagator as follows 




H 


(185) 


(186) 


where we denote 


Hf{x) 


tig 

E 


f{x + i) + f{x 
2ds 


*) 


(187) 


the operator H generates the hops for a random walk on a square lattice in dg 
dimensions. Using eq. (186), we can write the inverse propagator as 


2ds 1 




2d, <^0 




3gg ° a2 
2d, ^0 


1 ) 


-1 


+(^^0 + + i)-i 


H- 

(188) 


In this form, the inverse propagator appears as a sum of terms, each of which 

is an alternating product of H (i.e. single hops on the lattice that represents 

2 

space) and of the inverse of the operator dg + 1. We shall now rewrite this 
object in a way that clarihes its physical meaning. Let us start with its heat 
kernel representation 


Soo“ 


2d, "0 

The second step is the following identity 


—— = f dt exp (-1 (i + ) . 

d^ + l Jo ^ ^ ' ' 


— 8 

e2^ = 


r+oo 


dz 




e 2 a , 


or more explicitly 




r+oo 


dz 


V2' 


TTa 


e 2 a /(x + z) . 


(189) 


(190) 


(191) 
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In words, the operator exp(^i9^) is a diffusion operator that smears the target 
function by convolution with a Gaussian of variance a. But this interpretation 
is only possible if a is a positive real number. The closest formula if a is negative 
would be 

a 2 r+oo Ay 

fix) = p-| e 2|a| /(a; + i 2 ) , (192) 

but it requires that we complexify the variable x. This is precisely the situation 
we face when employing this formula to rewrite the inverse propagator: 
must be negative in order to obtain a properly normalized Gaussian. If we start 
from the Minkowski metric = +1), we can still get a Gaussian with the 
“correct sign” if we complexify the time. The conclusion of this digression is 
that the diffusive interpretation of the time evolution of a quantum system is 
only possible with imaginary time. 

From now on, we assume that = — 1. Eq. (189) can be rewritten as 


1 _ 

^ 2d, "o 


fixo) = 


dt 


d« 


r+oo 


27rta^ 


dz e * e fx-it /(xq + z) . (193) 


In the right hand side, we recognize an integral over diffusion processes of “du¬ 
ration” t, weighted by a factor exp(—1/2) that suppresses the contribution of 
diffusions longer than t ~ 1. Under such a diffusion, the time xq can shift by 
an amount of order a (the spatial lattice spacing). 

The interpretation of eq. (188) is now quite transparent. The inverse of 
is obtained by intertwining continuous diffusions in imaginary time (of arbitrary 
lengths, but weighted in such a way that the displacement in time is of order 
a) and single hops in one of the spatial directions. Note that in eq. (193), the 
integration over the length t can be done analytically by using 

rl+ 

/ 4 e-‘. (194) 

Jo yt 


Therefore, we can write 


1 - 



dz e f{xQ + z) , 


(195) 


where we denote 7 = . We see that the application of this operator is 

a convolution with a Laplace distribution. The typical width of the resulting 
smearing is z ~ 7“^ ~ a. In other words, the spatial lattice spacing also controls 
the typical size of the jumps in time that are interspersed between the jumps on 
the spatial lattice. This is consistent with our conclusion of the section 4, that 
the results become independent of the smallest lattice spacing when it becomes 
infinitesimally small (the continuous-time description adopted in this appendix 
corresponds to a discrete description of time with a time interval that goes to 
zero). In particular, the ultraviolet behavior in this limit is controlled by the 
next-to-smallest lattice spacing. 
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In the series expansion of the inverse propagator, we need to evaluate the 
n + 1-th power of this operator^®, acting on a starting distribution of the form 
(5(a:o), 


(^o) 


n+1 



(5(a:o + zi + ■ ■ ■ + Zn+i) 
1 n-l-l 


(196) 


For xq = 0, the evaluation of the integral is straightforward and one obtains 


a„(0) 


7 (2n)! 

2 4"n!2 ■ 


(197) 


The second factor is the probability that a random walk in one dimension returns 
to the origin after 2n steps. 

Using these results, we obtain the following expression for the leading term 
of(C(0)<^a(0)), 

(0* (0)<^(0)) = i^Co tradj (1) + • • • , (198) 


with 

Co = ^ y _^_ y __ ( 199 ) 

In this summation, we must have an even number 2n of spatial hops since their 
combination must form a closed loop on the spatial lattice. The first factor 
is therefore eq. (197) evaluated for 2n instead of n. The second factor counts 
the number of closed random walks of length 2n in dg dimensions. In order to 
disentangle the sums over ni, • • • , Ud,, we need a modified form of the Borel 
transformation, tuned to cancel the n-dependent prefactor. Notice Hrst that 


y r y ^ (yi 

V^Jo y® ^ 42" (2n)!- 


(200) 


This leads immediately to 


-1 j_ y 
2 y io y 



( 201 ) 


^®The choice of the Aq = 0 gauge is crucial here. Indeed, if we had a non-zero position 
dependent Aq background field, the hops in the time direction would not commute with the 
spatial hops, and it would be impossible to collect them as the n + 1-th power of a single 
temporal hop operator. 
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and finally to 


(r(0)</>(0)) 


_ tiadj (1) 1 f°° dr 

ad<oi. d-i-s-o -y/27rds Jq \/t 


rd. 




( 202 ) 


which is identical to the formula of the footnote 12, that was obtained by taking 
the limit —)■ 0 in a discrete time formulation. 
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